Temperature dependence of nonlinear auto-oscillator linewidths: Application to 

spin-torque nano-oscillators 
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The temperature dependence of the generation linewidth for an auto-oscillator with a nonlinear 
frequency shift is calculated. It is shown that the frequency nonlinearity creates a finite correlation 
time, r = l/r p , for the phase fluctuations. In the low-temperature limit in which the spectral 
linewidth is smaller than F p , the line shape is approximately Lorentzian and the linewidth is linear 
in temperature. In the opposite high-temperature limit in which the linewidth is larger than F p , the 
nonlinearity leads to an apparent "inhomogeneous broadening" of the line, which becomes Gaussian 
in shape and has a square-root dependence on temperature. The results are illustrated for the 
spin-torque nano-oscillator. 
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Spin-torque auto-oscillators (STOs), which are based 
on the spin-transfer effect in magnetic multilayered 
nanostructuresj^ possess a much stronger intrinsic fre- 
quency nonlinearity compared to conventional auto- 
oscillatory systemsi^ For example, variations in the gen- 
eration frequency can be two orders of magnitude greater 
than the generation linewidth in STOs, while this ratio 
is typically of the order of unity for conventional auto- 
oscillators. One important and interesting consequence 
of this nonlinearity is the strong coupling between the 
phase and amplitude fluctuations in the presence of ther- 
mal noise. While the amplitude fluctuations are subject 
only to a Gaussian white noise forcing, they act, through 
the nonlinear coupling, as a colored-noise source for the 
phase fluctuations^^ This coupling has been shown to 
cause substantial broadening of the generation linewidth 
above the generation threshold^ and to lead to asymmet- 
ric line shapes near threshold^ 

In this article, we investigate the temperature depen- 
dence of the generation linewidth of a nonlinear oscil- 
lator using a universal model of an auto-oscillator with 
a strong nonlinear frequency shift We demonstrate 
that the frequency nonlinearity creates a finite correla- 
tion time, r = l/r p , for the phase fluctuations caused 
by the thermal noise, which leads to a non-Lorcntzian 
shape of the power spectrum in general. A Lorentzian 
line shape is recovered only in the limit of low tempera- 
tures (at which linewidth magnitudes Au) are small com- 
pared to the inverse correlation time of the phase fluc- 
tuations Aw <C r p ). The generation linewidth in this 
regime is linearly proportional to temperature. In the 
opposite limit of high temperatures (at which linewidth 
magnitudes are large compared to the inverse correlation 
time Aw 3> r p ), the frequency nonlinearity leads to the 
apparent "inhomogeneous broadening" of the linewidth 
in which fluctuations in the oscillation amplitude leads 
to variations of the generated frequency by virtue of the 
strong frequency nonlinearity. The line shape in this limit 



is Gaussian with a linewidth that is proportional to the 
square-root of temperature. The results are illustrated 
for a spin-torque auto-oscillator. 

We study the stochastic dynamics within the frame- 
work of a universal nonlinear oscillator model with a 
stochastic force /„ that represents the action of thermal 
fluctuations) 5 ^ 7 ' 8 ' 9 

dc 

— + iu(p)c + r + ( P )c - r_( P ) c = f n (t), (l) 

where c = |c| exp(i<^>) is a dimensionless complex variable 
that describes the excited mode amplitude \c\ and phase 
4>, and p — | c | 2 is the dimensionless oscillation power. 
The equation of motion contains three terms which are 
functions of power: the nonlinear generation frequency 
Lu(p), the positive nonlinear damping T + (p) determined 
by the natural dissipative processes in the system, and 
the negative nonlinear damping T_ (p) determined by the 
active element that supplies external energy in the sys- 
tem. The stochastic term f n (t) describes the action of 
thermal fluctuations on the auto-oscillator and its statis- 
tical properties are determined from the thermodynamic 
properties of the oscillator at thermal equilibrium. This 
is achieved by taking /„ (t) to be a white Gaussian noise 
with zero mean, (f n (t)) — 0, with a second-order corre- 
lation function, (/ n = ^D n 5{t - t'), where the 
diffusion constant for a nonlinear oscillator is given by 7 - 



D n (p) =T + (p)r)(p) = r+(p) 



k B T 
\uj(p) ' 



(2) 



Here 77 is the effective noise power in the nonlinear regime 
and A is a coefficient that relates the oscillator power p 
to the oscillator energy £. 

In the following, we use the STO as an example of 
the stochastic nonlinear auto-oscillator described by the 
above model. At reasonably small generation powers, 
p < 1, the coefficients of the power-dependent terms in 
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Eq. JTJ) can be approximated with the following Taylor 
series expansions^ w (p) = uj +Np, T + (p) — T (l+Qp), 
and r_(p) = <rl(l — p), where I is the bias current and 
the explicit definitions of the coefficients N, Q,<J, Tq are 
given in Ref. [y. The coefficient of proportionality be- 
tween the energy and power of the STO is A = V^ffMo/7, 
where V c s is the effective volume of the magnetic ma- 
terial of the free layer involved in the auto-oscillations, 
Mo is the saturation magnetization of the free layer, and 
7 is the gyromagnetic ratio. The stationary generation 
power of the STO is given by p = (( — !)/£ + Q), where 
( = I/Ith is the supercriticality parameter. The thresh- 
old current for auto-oscillations, Ith = Fo/er, is deter- 
mined by the condition by which the linear components 
of the negative and positive damping terms are equal, 
F_(0) = r + (0). In the supercritical region, sufficiently 
far above the generation threshold, the stationary power 
Po of a nonlinear auto-oscillator is found from the condi- 
tion r + (p ) = T-(p ). 

To determine the generation linewidth in the far su- 
percritical regime, it is sufficient to linearize Eq. ([1]) by 
examining the time evolution of the power fluctuations 
6p(t) = p(t) — po about the steady-state trajectory, as Sp 
is much smaller than the stationary value po. By lineariz- 
ing Eq. JT]) near the point p = po, we find the following 
two equations for the power dp and phase 4> fluctuations 
of the auto-oscillations, 



correlation function of the form, 



dSp 
~~dt 



2Y v 5p = 2v5£Re[/„(t)], 



■u)(p 



Im[f n {t)]-NSp, 



(3a) 
(3b) 



where f n (t) = f n {t)e~ l ^ tS> is a stochastic process with the 
same properties as f n (t), N = d p uj(p) is the coefficient of 
the nonlinear frequency shift, 



T p = d p (T + - r_)p = G eB po, 



(4) 



is the damping rate for the power fluctuations, and all the 
derivatives are calculated at p = pa . It is immediately ev- 
ident that the nonlinearity N leads to a coupling between 
the phase and amplitude fluctuations, which generates 
an additional noise source —NSp for the phase variable 
(j). The physical mechanism behind this additional noise 
term can be understood in terms of an additional fre- 
quency modulation uj{p{t)) ~ w(po) + N Sp which arises 
from the power fluctuations Sp(t) and the frequency non- 
linearity. This additional noise term can be considered as 
an inhomogeneous broadening of the oscillator linewidth: 
oscillators with different powers (due to thermal fluctua- 
tions) p have different oscillation frequencies u>(p). 

Since the stochastic system Q is a linear system of 
equations and the noise f n {t) is a Gaussian process, both 
Sp(t) and <j>{t) are also Gaussian processes for which the 
complete set of statistical characteristics can be easily 
obtained. As such, we find a vanishing mean value for 
the power fluctuations, (5p(t)) — 0, with the two-time 



(Sp(t)Sp(t')) =p v(Po 



r+(p ) 
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(5) 



The mean value of the phase <j)(t) (for a given initial value 
ipQ at t = 0) is given by (4>{t)) = (f>o — u>{po)t, and the 
variance of the phase fluctuations is given by 



A(j) 2 (t) = 2Au 
where 



(l + v 2 )\t\-v- 



1 _ e -2r p |t| 



2r c 



dpW 



N 



d p (T+ - r_) G c ff 



(6) 



(7) 



is the normalized dimensionless nonlinear frequency shift 
coefficient which characterizes the relative influence of 
the frequency nonlinearity N compared to the nonlinear- 
ity of total effective dissipation G e g on the auto-oscillator 
behavior, and 2Aujq is the generation linewidth of a linear 
(y = 0) auto-oscillator, 



2Aw = 



k B T 



(8) 



In stark contrast to the case of a linear oscillator for 
which \i/\ <C I, the phase variance A(j) 2 (t) for an auto- 
oscillator with strong frequency nonlinearity v 3> I has 
a nonlinear dependence on the time interval \t\. This 
nonlinear dependence results from the finite correlation 
time r = l/r p of power fluctuations and by the addi- 
tional nonlinear source of colored noise —Ndp(t) in the 
equation of motion for the phase variable ([3]) as men- 
tioned previously. Only for time intervals \t\ that are 
much larger than the correlation time of the power fluc- 
tuations I /r p does the phase variance vary linearly as a 
function of time. Thus, the power spectrum of a nonlin- 
ear auto-oscillator is non-Lorentzian in general. 

The power spectrum S(Q) of the nonlinear auto- 
oscillator ([T]) in the presence of thermal noise, 



s{n) = j ' dr jc(t) ( 



(9) 



is, up to a constant of proportionality, the Fourier trans- 
form of the two-time correlation function of the oscillator 
variable c(t), 



IC(r) = (c(t + r)c*(t)). 



(10) 



In the supercritical regime in which the oscillator power 
is stable and phase fluctuations dominate, we can con- 
sider simply the phase noise contributions to the power 
spectrum. In this case, the autocorrelation function cor- 
responding to the random phase can be written as 

JC{t) = Po e Mt)-+m exp [-A0 2 (t)/2] . (f f ) 

For an STO, the power spectrum measured in a typi- 
cal experiment corresponds to the Fourier transform of 
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the autocorrelation function of the variable voltage that 
arises from variations in the giant or tunnel magnetore- 
sistance. The relation between the autocorrelation func- 
tion of the STO voltage and the autocorrelation function 
JC(t) of the spin- wave auto-oscillation has been derived 
elsewhere (e.g., see Eq. (10) in Ref. 0). 

The form of the phase variance in ^ does not permit 
a simple analytical expression for the Fourier transform 
of the autocorrelation function (fTTj) to be obtained in 
the general case. However, if the temperature T and, re- 
spectively, generation linewidth Aoj are sufficiently small, 
Aoj <C r p , one can neglect the exponential factor in © 
at the characteristic decoherence time scale t ~ 1/Aw, 



where the characteristic linewidth is defined by 



A<t> 2 (t) w 2Aw (l + v 2 )\t\ 



v 

2TV 



(12) 



In this low-temperature limit the autocorrelation func- 
tion IC(t) is a simple exponential function, which leads to 
a power spectrum <S(f2) with a Lorentzian shape and a 
full width at half maximum of 5 



2Acc>lt = 2Acj (1 



r+bo) 



k B T 
£(po) 



(1 + ^ 2 ). (13) 



In this limit, the frequency nonlinearity leads to an in- 
crease in the generation linewidth by a factor of (1 + v 2 ) 
with a linear dependence on the temperature. One im- 
mediately observes that the condition Aw <C F p used to 
derive Eq. <|13|) can be rewritten as a condition for the 
temperature 



k B T<^ 



r+(po) 



£(po) 



(14) 



or expressed as a condition for the relative magnitude of 
the power fluctuations, Ap/po <C (1 + v 2 )^ 1 ! 2 . Thus, 
(|13p represents the low-temperature asymptotics of the 
generation linewidth for a nonlinear auto-oscillator. Es- 
timates made using Eq. (|13p for typical parameters of 
spin-torque oscillators (permalloy circular nanopillar of 
the radius R c = 50 nm and thickness of the free layer 
L = 5 nm, generation frequency w(po)/27r ~ 30 GHz) 
show that (|13p is quantitatively correct for temperatures 
in the range T < 10 — 100 K, depending on the super- 
criticality £ (see, e.g., Figs. 2a and 2b in Ref. @). 

A simple analytical solution for the linewidth can also 
be obtained for the opposite limit of high temperatures 
in which the correlation time of the power fluctuations 
(l/r p ) is much longer than the coherence time of the os- 
cillator (~ 1/Aw). In this limit, the exponential function 
in ([6]) can be expanded in a Taylor series to give 



Acj> 2 (t) « 2Awo(|f| + v 2 T p t 2 ). 



(15) 



For large frequency nonlincarities, v 3> 1, one can neglect 
the linear term in \t\ and keep only the quadratic term 
in the Taylor series expansion, which leads to a Gaussian 
power spectrum, 



S(fi) ~ exp 



(w(po) - n) 2 



2Aw| T 



(16) 



AciJht 



r+(p )r F 



I k B T 

£(poY 



(17) 



It is interesting to note that the generation linewidth 
in this regime is proportional to yT, in contrast to the 
linear temperature dependence (| 13[) of the linewidth in 
the low-temperature regime. Such a temperature depen- 
dence has been originally proposed in Ref. 03 based on 
macrospin simulations of an STO. In the context of the 
present work, this linewidth behavior can be understood 
as an apparent inhomogeneous broadening by power fluc- 
tuations which are translated into a modulation of the 
oscillation frequency through the frequency nonlinearity. 
The expression (II T[) for the generation linewidth is valid 
in the temperature range 



r +{Po 



£(po) 



r +{Po 



£(po), (18) 



which can also be rewritten in terms of the range for the 
power fluctuations as \/\v\ <C Ap/p <C 1. 

To study the transition between the low and high tem- 
perature asymptotic behavior, we performed numerical 
calculations of the Fourier transform of autocorrelation 
function (fTTj) using the phase variance given by the gen- 
eral expression For the numerical work it is conve- 
nient to use the dimcnsionlcss definitions of time, t' — 
r p t, and temperature, T = [(T + ( Po )k B )/(T p S (p ))] T. 
With these scaled variables, the phase variance in ([6]) 
can be simplified to 

A<j> 2 (t') = T' \(1 + v 2 )\t'\ - (v 2 /2){l - e 21 *' 1 )] . (19) 

The results of the numerical calculations of the genera- 
tion linewidth (normalized by F p ), for several values of 
the normalized nonlinearity coefficient is, as a function of 
the dimensionless temperature T 1 over 8 orders of mag- 
nitude are shown in Fig. [1] The linewidth is determined 
from the full width at half maximum of the power spec- 
trum, without fitting to any particular line shape. 

In the absence of nonlinearity the linewidth is sim- 
ply proportional to the temperature, and in our reduced 
units, there is a one-to-one correspondence between the 
two quantities. As the nonlinearity v is progressively in- 
creased, the temperature range over which the linear be- 
havior holds decreases over the temperature range stud- 
ied, or, in other words, the linear regime is pushed to- 
ward lower temperatures, as indicated by (|14p. In the 
low temperature regime the linewidth dependence on v 
is consistent with Eq. (fl"3"]) , i.e. the linewidth increases as 
(1 + is 2 ), as expected. As the temperature is increased, a 
gradual transition toward the square-root temperature 
dependence is observed. Over the temperature range 
studied, this transition occurs more rapidly for the larger 
values of v, as expected from (fT8"|) . The low and high tem- 
perature limits are shown explicitly in Figs. [T}) and [TJ;, 
respectively. 
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FIG. 1: (Color online) Normalized linewidth, Auj' = Aw/r p , 
as a function of the dimensionless temperature T' . (a) 
Linewidth variations for different values of the frequency non- 
linearity v. (b) Low temperature asymptotics for v 2 — 100 
(circles are numerical data, solid line is fit using the linear 
function A/' = 100T'). (c) High temperature asymptotics 
for v = 10 4 [circles are numerical data, solid line is a fit using 
the function A/' = 222 (T') 518 ]. 



It is possible to test our theory directly in experiment 
by tuning the nonlinearity coefficient v over a reason- 
ably wide interval. For STOs, this can be achieved by 
varying the orientation of the out-of-plane external bias 
magnetic field relative to the plane of the "free" layer for 
magnetic nano-contact system^ or by varying the orien- 



tation of the in-plane bias magnetic field relative to the 
"easy" anisotropy axis in an anisotropic in-plane magne- 
tized magnetic nanopillarJ^ In both cases, the variation 
in v is due to different dynamic dipolar fields generated 
as a result of changes in the magnetization precession 
axis<2 We would like to stress that any quantitative com- 
parison between our theory and experiment also requires 
detailed knowledge of the temperature dependence of all 
the key magnetic parameters of the STO, such as oscil- 
lation frequency, shape anisotropy, and relaxation rate, 
so that the supercriticality can be determined accurately 
for a given value of the bias current. One way to over- 
come such difficulties is to drive the STO with an external 
noisy microwave magnetic field at fixed temperature. In 
this case the effective noise temperature is proportional 
to the power P n of microwave noise source and the de- 
pendence Auj(P n ) should be qualitatively the same as the 
temperature dependence Au(T), in particular, Auj ~ P n 
for small P n and Auj ~ \[P^ for relatively large P n . 

In summary, we have developed a theory for the tem- 
perature dependence of a generation linewidth of an auto- 
oscillator with a large nonlinear frequency shift. At low 
temperatures, the nonlinearity leads to the renormaliza- 
tion of the phase noise, which results in a Lorentzian 
power spectrum and a linear temperature dependence. 
At high temperatures, the nonlinearity leads to an ap- 
parent inhomogeneous broadening of the spectral line, 
giving rise to a Gaussian line shape and a square-root 
dependence on the temperature. 
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